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We investigate the low-temperature electrical and thermal transport properties in atomically pre- 
cise metallic heterostructures involving strongly-correlated electron systems. The model of the Mott- 
insulator/ band- insulator superlattice was discussed in the framework of the slave-boson mean-field 
approximation and transport quantities were derived by use of the Boltzmann transport equation in 
the relaxation-time approximation. The results for the optical conductivity are in good agreement 
with recently published experimental data on (LaTiOa)jv/(SrTi03)M superlattices and allow us to 
estimate the values of key parameters of the model. Furthermore, predictions for the thermoelectric 
response were made and the dependence of the Seebeck coefficient on model parameters was studied 
in detail. The width of the Mott-insulating material was identified as the most relevant parameter, 
in particular, this parameter provides a way to optimize the thermoelectric power factor at low 
temperatures. 

PACS numbers: 71.10.Fd, 71.27. +a, 72.10.-d, 73.21.Cd 



I. INTRODUCTION 

Artificial heterostructures and superlattices composed 
of different perovskite oxides have received a considerable 
attention in recent years. In many cases the heterostruc- 
tures are based on two insulators that, interestingly, ex- 
hibit metallic behavior in atomically precise superlat- 
tices. For instance, Ohtomo and coworkers reported 
on a metallic conductivity in LaTi0 3 /SrTi0 3 (Ref. [l|) 
and in LaA10 3 /SrTi0 3 (Ref. [3) heterostructures. Thiel 
et al£ demonstrated the possibility to tune the carrier 
density in the latter system by electric-ficld-effect, and 
Reyren et alA detected a superconducting transition at 
T c = 200 mK. More recent studies involve LaV0 3 /SrV0 3 
heterostructures^ and LaV0 3 /SrTi0 3 interfaces^ Simi- 
lar to the LaTi0 3 /SrTi0 3 system, the latter system falls 
into the general class of band- insulator/ Mott-insulator 
heterostructures. 

These experimental findings have stimulated a consid- 
erable amount of theoretical workjl~— Common to all 
these studies is the technical challenge to handle strong 
local correlations in a spatially non-uniform system and 
therefore calculations are mainly based on effective mod- 
els. Similar to interfaces between semiconductors, elec- 
tronic charges are redistributed in order to maintain an 
electrostatic stable solution. In the case of two insulating 
materials, the mutual doping can lead to metallic behav- 
ior of the interface. Generally, its electronic phase may 
differ from the bulk phases of the two constituents, and 
a variety of ordered phases in charge, spin, and orbital 
degrees of freedom, as well as more exotic phases, are 
predicted. 

In this work we are exclusively concerned with the ex- 
perimentally observed metallic properties ^ 19 ' 20 of the 
interface between Mott insulator (MI) and band insula- 
tor (BI). From the theoretical point of view, there has 
been relatively little attempts to describe the transport 
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FIG. 1: (Color online) Schematic view of the setup considered 
in this article. N is the number of unit cells of the Mott insu- 
lating material (MI) and M is the number of unit cells of the 
band insulating material (BI). In addition, the conduction- 
band density is indicated. We assume a thermal gradient VT 
or an electrochemical field E along the x direction. 



properties of these quasi-two-dimensional metallic elec- 
tron systems and to clarify the role of the expected strong 
electron-electron interaction. It is therefore of interest to 
study some general aspects of the electronic transport on 
a qualitative level. We show that experimental data2°- 
on the optical conductivity of LaTi0 3 /SrTi0 3 superlat- 
tices were well reproduced for reasonable parameters of 
the present theory. In addition, we focus, in this work, 
on the low-temperature thermoelectricity 

A schematic view of the considered setup is shown in 
Fig. [1] We assume a cubic lattice structure with per- 
fect lattice match and a superlattice modulation in the 
z-direction. This setup is motivated by experiments on 
(LaTi0 3 )iv/(SrTi0 3 )j\f superlattices but other material 
combinations may be possible. Quite in general we as- 
sume a polar perovskite Mott insulator A'B0 3 with a 3d 1 
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configuration and a nonpolar perovskitc band insulator 
ABO3 with a 3d configuration. The transport prop- 
erties were calculated in response to an electrochemical 
field E or a temperature gradient VT in the x direction. 
This is in contrast to previous theoretical work where 
the transport properties in strongly-correlated multilayer 
nanostructurcs along the z direction was studied*^ 

The motivation to study the thermoelectricity is based 
on several interesting observations: (i) The solid solution 
Sr 1 _ 2 ,La 2; Ti03 shows a large thermoelectric response 2 ^ 
for < x < 0.1. (ii) From the study of narrow-band or- 
ganic conductors^ intermetallic compounds^ - — as well 
as cobalt or other transition-metal oxides^ it is known 
that correlation effects can enhance thermoelectricity, 
(iii) Lower-dimensional structures, such as quantum-well 
supcrlatticcs, offer additional parameters to optimize the 
thermoelectric response and, at best, overtop the bulk 
properties. 27 ' 28 

This paper is organized as follows: In Sec. II we in- 
troduce the model. In Sec. Ill we discuss a quasiparticlc 
approach to the low-energy properties and in Sec. IV the 
transport coefficients are calculated. A short overview of 
the obtained results is given in Sec. V with a subsequent 
discussion in Sec. VI. We summarize our main conclu- 
sions in Sec. VII. 



electron and electron-counterion interactions as follows^ 
H = -t^2 + u 5Z"*T"a 

(ij)a i 

+ Vfri + \Y1 rnWijhj. (1) 

The countcrions sit in between the electronic sites, and 
simulate the electrostatic difference between the A (Sr 2+ ) 
and the A' (La 3+ ) cations. The electrons at site i feel an 
electrostatic potential 

*=-^£jf^=r (2) 

where r^ on denotes the position of the counterions. The 
interaction matrix was given by Wij = Ec /\f*% — Tj |, and 
we focus on charge neutral systems. We neglect effects of 
the lattice relaxation and fix the counterions at equidis- 
tant positions. As discussed in Ref. [H it was expected 
that the screening provided by the atomic reconstruction 
reduces the conduction-electron density on the central 
layers and enhances the density in the band insulator 
away from the interface. The modeling of such effects is, 
however, beyond the scope of this paper. 



III. QUASIPARTICLE DESCRIPTION 



II. MODEL 



We adopt a single-band model, thereby neglecting 
complicating aspects regarding the orbital degrees of free- 
dom typically present in transition metal oxidesi 2 ^ In this 
model three different energy scales are important: (i) the 
hopping matrix element t, (ii) the onsite Hubbard U, and 
(iii) the Coulomb parameter of the long-range interac- 
tion Ec = e 2 /££>a, involving the dielectric constant Ed 
of the core electron background and the lattice constant 
a w 3.9AJ£ The onsite Hubbard U is intrinsically large, 
namely U > U c where U c ~ Uf R is the critical interac- 
tion strength for the bandwidth controlled Mott transi- 
tion in the bulk system and U R R « 16t is the Brinkman- 
Rice valued However, it is instructive^ 7 - at some points 
in the following discussion to vary the value of U from 
the so-called Hartree regime U <C U c to the Mott regime 
U > U c in order to see how the electronic properties 
changed A realistic value of the Coulomb parameter 
is Ec = 0.8t. A rough estimat e of t he Thomas-Fermi 
screening length gives Xtf ~ \ftjE~c~a, which is of the 
order of the lattice constant a. In the limit Ma 3> Xtf 
it is, in a first approximation, sufficient to consider only 
one quantum well made out of N layers of the MI ma- 
terial embedded in two semi-infinite BI systems. The 
microscopic model for the quantum-well system is de- 
fined by an extended single-band Hubbard Hamiltonian 
on a simple cubic lattice including long-range clcctron- 



A. Single-particle Green's function 

In the following we consider the low-temperature limit 
T — > 0. We apply a quasiparticlc description for the low- 
energy properties of the above model system. Follow- 
ing Refs. I32I andl33l the layer-dependent Green's function 
reads 



G KV (k,w) = [u + m - t(k) - ^(k.w)]^, 1 



(3) 



where the matrix i(k) is the Fourier transformed hopping 
matrix 



t(k) = 




(4) 



and I, I' subscript the layers. The free dispersion of the 
two-dimensional nearest-neighbor hopping tight-binding 
model is £k = — 2£(cos k x a + cos k y a) and S ff (k, u>) is the 
layer-dependent self-energy. We assume that the self- 
energy is local, thus independent of k, and diagonal in 
the layer index. Furthermore, for ui — > we assume a 
paramagnetic Fermi-liquid form 



Sh'(cj) = 5 W 



A, 



(1 



1 



(5) 



where we have neglected the imaginary part, which is 
expected to be proportional to u) 2 and have suppressed 
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the spin index a. From Eq. ((5|) it follows that the layer- 
dependent renormalization factor, 



is purely local, 



\ dt{io) 










u)=0- 


Z w = du'Z 
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(6) 



(7) 



A; in Eq. ([5]) is an effective layer-dependent chemical po- 
tential including Hartree-Fock like (electrostatic) contri- 
butions from the long-range Coulomb potential as well as 
lu = contributions from the dynamical self-energy due 
to the onsite repulsion. The low-energy part of Green's 
function may be written in terms of an effective quasi- 
particlc Hamiltonian 



G(k,w) = VZ ^-# efr (k)J VZ + G inc (k,u>), (8) 

where the effective Hamiltonian is given by 

^i^k + Ai — tziz 2 
-tz\Zi zfe k + A 2 -tz 2 z 3 
effW_ -tz 2 z 3 z|e k + A 3 ~tz 3Zi 



and Gi„ c (k, lu) denotes the incoherent part of Green's 
function, which we will neglect for lu — > 0. Let us intro- 
duce the eigenstates |kf) of if e fj(k) as 



fleff (k)|M = E kv \kv), 



(9) 



where v is the subband index and the envelope function 
is given by ip^{l) = (Z|ki/). The k dependence enters via 



Eventually, we find for lu — > 



G(k, w ) 



w — E kv 

Green's function for the quasiparticlcs is 



(10) 



(11) 



G<2 p (k, 



E 



and is diagonal in the quasiparticle subband basis. 
B. Slave-boson mean- field approximation 



(12) 



suggested in Ref. [l7| the use of a layer-dependent gen- 
eralization of the Kotliar-Ruckenstein slave-boson mean- 
field (SBMF) approximation.— In the following, when it 
comes to actual numerical results, we rely on this ap- 
proximation for the low-temperature electronic proper- 
ties. (For an extended discussion on the validity of this 
approach see, e.g., Refs. and HE) The SBMF ap- 
proximation was successfully applied in other contexts to 
study spatially non-uniform solution a 37 ! 38 and compares 
well with (two-site) DMFT calculations in the present 
context . 15 ! 17 In the SBMF approximation, one directly 
addresses the properties of the quasiparticles by con- 
structing a free energy, which takes into account the 
competition between kinetic and potential energies. It 
also allows to account for local electron-correlation ef- 
fects on the low-energy properties caused by the strong 
onsite repulsion. In particular, z\ and A/ are determined 
self-consistently in a numerically efficient way. 

The results of the SBMF approximation to the 
quantum-well model [Eq. ([I}] were described in detail in 
Ref. [I?] and we will recapitulate only some basic aspects 
needed for a general understanding of the present work. 
Within this approximation the long-range Coulomb in- 
teraction is treated in a Hartree-like way and the onsite 
interaction in the spirit of the Gutzwiller approximation 
as known from the homogeneous Hubbard model (see, 
e.g., Ref. l39l) . We introduce the parameter U r = U — Ec- 
This parameter is chosen such that for U r = the 
Hartrcc mean-field solution is recovered. We then deter- 
mine self-consistently the electron-density profile in the 
z-direction, the screened effective potential A;, and the 
fraction of doubly occupied sites in each layer I. The 
self-consistent treatment allows us to find the quasipar- 
ticle dispersion E ku and the envelope function ipicv(l) of 
the quasiparticle |ki/). Due to the strong onsite Coulomb 
repulsion, the quasiparticle velocity is renormalizcd, 



Vki/ = —z- V£k- 
n 



The renormalization amplitude Z kl/ is given by 



(13) 



de k 



< 1. 



(14) 



It depends both on v and £k- In particular, the Fermi 
velocity of the subband v is reduced by 



Z v 



dE u 



ds-k 



(15) 



The low-energy behavior of the self-energy may be ob- 
tained using the dynamical-mean-field-theory (DMFT) 
approximatio n 7 i 14 i 15 from which the effective quasipar- 
ticle Hamiltonian ff e g(k) can be obtained. As a nu- 
merically less expensive alternative the present authors 



For quasiparticles whose wave functions involve large 
contributions from the localized orbitals in the center of 
the Mott-insulating material, Z v is considerably smaller 
than one, meaning that near the chemical potential the 
resulting subbands are almost flat. 
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IV. CALCULATION 
A. Transport coefficient 

The self-consistent renormalizcd quasiparticle sub- 
bands at T = were taken to calculate the transport 
properties of the two-dimensional electron system. This 
restricts our analysis to lowest order in T. 

We adapt the two-dimensional version of the Boltz- 
mann transport theory and consider the effect of the 
electron-electron interactions as reflected in the renor- 
malization of the quasiparticle velocity [Eq. (fT4|) ] and in 
the k dependence of the envelope function From 
a Fermi-liquid viewpoint^ this means that we neglect 
the (residual) interactions between Landau's Fermi-liquid 
quasiparticles. The relation to the linear-response Kubo 
formula is discussed in Appendix [SJ We introduce the 
local quasiparticle distribution function /k„. The lin- 
ear response to an applied uniform electric (electrochem- 
ical) field E(w) or temperature gradient VT in the di- 
rection perpendicular to the growth direction of the het- 
erostructure (xy-plane) is found by linearizing the dis- 
tribution function /k„ = /°(£'ki/) + .9ky, where f°(E) = 
[1 + exp(/3i?)] _1 is the equilibrium Fermi-Dirac distribu- 
tion function for the inverse temperature (3 = 1/kgT and 
where gu v is proportional to the applied field. 

For the setup shown in Fig. [T] it is sufficient to restrict 
the calculations to the x component of the currents. In 
this case we define the transport distribution function^* 
<&{E) as follows: 



<Z>(E) = 



Tu(E) 



iijjT v {E) 



v v {E)S v {E). (16) 



Here, we introduced the relaxation time t v (E) which is 
assumed to be independent of the two-dimensional mo- 
mentum k. S V {E) is the area of the constant energy 
surface at energy E of the subband v and 



v v {E) 



1 



S V {E) 



\v kv \dS y (E) 



(17) 



is the averaged velocity over the constant energy surface. 
According to Eq. (Tl6|) the transport distribution func- 
tion is a sum of the contributions of each partially filled 
subband, 



(18) 



following Fermi integrals over the single kernel function 
$(£), 



(20) 



For the conductivity o- 2 Di the thcrmopower or Seebeck- 
cocfficicnt Q2D1 an d the thermal conductivity ftf^, one 
then find a 41 ' 42 



o w = L<°>, 



i2D = 



_1_L« 

'eTL(°)' 



and 



V 2D 



e 2 T 



L (2)_ 



(21) 
(22) 

(23) 



Note that the thcrmopower Q does not contain a volume 
term and therefore we approximate the thermopower of 
the superlattice as 



Q-Q21 



(24) 



On the other hand, electrical and thermal conductivities 
of the superlattice arc approximated by 



0~2D 



a(M + N) 



and 



V 2D 



a(M + N) ' 



(25) 



B. Impurity scattering 



In the following we will estimate the relaxation 
time t v (E) entering the transport distribution function 
[Eq. HH)]. We assume that at low temperatures the dom- 
inant relaxation mechanism is clastic scattering by im- 
purities or vacancies and in the following we restrict the 
discussion to this case. For simplicity we consider the 
short-range impurity Hamiltonian 



(26) 



where i' labels the lattice sites of the impurities. We ne- 
glect multiple-scattering by different impurities and use 
the single-site T-matrix approximation 4 ^ with an aver- 
age over the locations i of the impurity. This average 
introduces a sum over layers I and we find 



Alternatively, we introduce the layer-resolved transport 
distribution function as 

= J2$v(E)i/>Ev(l) 2 , $(E)=J2M E )- (19) 



The notation ipEv(l) is short-hand for tp^Q) with E^ v = 
E. Equations (|18p and l|19p allow us to study the different 
contributions to the transport coefficients. We define the 



E 

1 



l-V G Q E P (l,l) 



(27) 



Here, T^J, is the atomic T-matrix for a single impurity 
at site i and 



(2tt) 2 E - Euu + «0+ 



(28) 



5 



is the retarded local quasiparticle Green's function of the 
pure system. We neglect the real part of G® P and intro- 
duce the layer-resolved quasiparticle density of states 



M = 10 



p(l,E) 



1. 

n 



— — Im G QP 



E 



E (hi) 
d 2 k 



(2nf 



(29) 



The relaxation time follows as 



Tu(E) 



h V 2 c imp Y / I 



1> El/ (l) 2 p(l,E) 

+ 7T 2 vi P (i,Ey 



(30) 



where 



•imp 



is the impurity concentration assumed to be 



independent of the layer I. 

We note that even though the energy dependence of 
the relaxation time near the Fermi energy is rather weak, 
varying the potential strength Vq has a significant in- 
fluence on the transport properties and, in particular, 
on the scaling behavior of these quantities with N (see 
Sec. El). The structure of Eq. §U§ offers the possibility 
to smoothly vary the relaxation time in each subband 
and to study different scenarios of quasiparticle scatter- 
ing. Due to the spatial non-uniformity, the value of p{l, 0) 
strongly changes with / reaching its maximum in the cen- 
ter of the heterostructure where the quasiparticlcs states 
are confined to energies near the Fermi energy E = 0. 
The first Born approximation is recovered in the limit of 
weak impurities: ttVop(1,0) -C 1 for all I. In this case, 
t„(0) is substantially smaller for the subbands whose spa- 
tial weight is mostly located in the center of the quantum 
well. In other words, slower quasiparticlcs arc scattered 
more strongly. In this case, the quasiparticle transport 
is mainly determined by the properties of the almost 
empty subbands. In the unitary limit, TrVop(l,0) 3> 1, 
the scattering rate of the slower quasiparticlcs was sig- 
nificantly reduced. In this case, the quasiparticle trans- 
port is mainly determined by the strongly renormalized 
subbands. For moderate values of Vq, the dependence of 
the scattering rate on the subband index was less pro- 
nounced. 

Beside the impurity potential strength, there arc other 
parameters that can influence the subband dependence 
of r„ in real systems. For example, it is conceivable that 
the scattering occurs mainly due to imperfections of the 
interface. Furthermore, there is no reason not to believe 
that at finite temperatures electron-electron or electron- 
phonon scattering will affect the individual subbands in 
a different way. In this sense, Eq. ([50)) can be considered 
to be a convenient form to vary the subband dependence 
by a single parameter. 




0.2 0.3 

N/(N + M) 

FIG. 2: (Color online) The plasma frequency uj p o of 
(MI)jv/(BI)m=io superlattices, as found by the present the- 
ory for t = 0.3 eV, and U r = 22t (dotted line), t = 0.15 eV 
and U r = (dashed line) and t = 0.15 eV and U r = 22t 
(solid line) . The filled circles denote the measurements of Seo 
et al. on (LaTiO3)jv/(SrTiO3)M=i0 superlattices reported in 
Ref . El 



V. OVERVIEW 
A. Optical conductivity and free carrier response 

Figure [5] shows a comparison of the calculated oscilla- 
tor strength lo p d of the Drude peak in the optical con- 
ductivity with the results of optical experiments reported 
by Seo et al. in Ref. HO on (LaTi0 3 ) JV /(SrTi0 3 )M=io 
superlattices. There is a good qualitative agreement be- 
tween the prediction from the theory in the Mott regime 
U r > Uf R « 16t and the measurements. Even a rather 
good quantitative agreement is achieved by choosing the 
hopping matrix clement of the order t = 0.15 eV. In con- 
trast, in the Hartree regime U r <C Uf the qualitative 
behavior of lo p d as function of N/(N + M) is quite dif- 
ferent and clearly not consistent with the experimental 
data. 

In Sec. IVII we also discuss the possibility to gain infor- 
mation about the scattering parameters by analyzing the 
width of the measured Drude peak. 



B. Thermoelectricity 

Due to the uncertainty in the subband dependence of 
the relaxation time t v it is difficult to estimate the ther- 
moelectric transport even at low temperatures. As an il- 
lustrative example we compare the results using Eq. (|30|) 
in the first Born approximation and in the unitary limit. 

We first show the results of the Born approximation. 
Figures 02 |4] and \5\ show the calculated low-temperature 
thermoelectric power (Seebeck coefficient) Q, the elec- 
tronic contribution K e to the thermal conductivity, and 
the power factor PF= oQ 2 as function of temperature 
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FIG. 3: The temperature dependence of the thermopower Q 
for quantum wells with N = 1,...,8, using the first Born 
approximation for the scattering process. The dotted-dashed 
line indicates the result for the N = 5 heterostructure when 
taking the self-consistent renormalized subbands at finite tem- 
peratures. The values of the relevant parameters are indicated 
in the figure. 
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FIG. 5: The power factor PF=aQ 2 for the same parameters 
as in Fig. [4] The dotted-dashed line has the same meaning as 
in Fig. |1 
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FIG. 6: The thermopower Q at k B T = O.Olt for N = 1 - 8 in 
the limit of a strong impurity potential Vo = Wt. 



FIG. 4: The electronic contribution K e to the total heat con- 
ductivity for the same parameters as in Fig. [3] The dotted- 
dashed line indicates the result for the N = 5 heterostructure 
when taking the self-consistent renormalized subbands at fi- 
nite temperatures. 



T for different width TV of the MI region, respectively. 
All the calculated transport quantities show the typical 
metallic low-temperature dependence. In particular, the 
thermopower Q is negative and its magnitude increases 
linearly with T. The numerical values of the relevant pa- 
rameters were estimated by comparing the theoretical re- 
sults to transport^ and spectroscopic^ measurements of 
SrTiC>3/LaTi03 superlattices and are indicated in Figs. [3] 
-[5] In order to check the validity of the zero-temperature 
approximation for the quasiparticle subbands, we also 
show the results for the N = 5 heterostructure as dotted- 
dashed line when taking the self-consistent renormalized 
subbands at finite temperatures. The difference between 
the two results is of order T 2 and in the considered tem- 
perature range, it is a good approximation to use the 



zero-temperature subbands for the calculations. 

In the unitary limit, a completely different behavior of 
the low-tempcraturc thermopower is observed. In Fig. [S] 
the Seebeck coefficient is shown for different width N of 
the quantum well. Interestingly, the value of the low- 
temperature thermopower Q shows an even-odd oscilla- 
tion with a decreasing amplitude for N — > oo. 

These results show that the Seebeck coefficient is in 
fact a sensitive probe of the scattering mechanism. Mea- 
surements of this quantity could reveal important infor- 
mations about which subbands contribute dominantly to 
the transport. We present a more detailed analysis on 
this point in the next section. 



VI. DISCUSSION 

A. Optical conductivity 

The optical conductivity is considered in the regime 
lot v ^> I of Eq. (|16[) . In this regime the conductivity 
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FIG. 7: (Color online) The averaged relaxation time f as 
found in the present theory (t = 0.15 eV) and U r = 22t for 
different scattering parameters as function of N. The filled 
circles denote the inverse width of the measured Drude peak in 
(LaTi03)jv/(SrTi03)M=io superlattices reported in Ref . I20I. 



In this expression D = D(N, U, Eq, t) is a function of the 
quantum well parameters. The calculations of lo p d using 
Eq. (|35|) arc shown in Fig. [2] In order to gain informa- 
tion about the scattering parameters, we assume that all 
the samples are of equal quality, Ci mp (iV) = const., and 
compare in Fig. [7] the value of f as a function of N with 
the width of the measured Drude peak. For a moderate 
impurity strength Vq ~ 1.8i and a reasonable impurity 
concentration c; mp « 0.085, the scaling of f with N, as 
well as its absolute value, is in good agreement with the 
experimental results. For these parameters the mean re- 
laxation time is f « 40 fs. 



B. DC conductivity 

The electrical and thermal conductivities were ob- 
tained in the static limit lot v <C 1 of Eq. (flU)) . From 
Eq. (|2"Tj) it follows that at low temperatures the dc con- 
ductivity is constant and given by 



o"2d(w) at low temperatures has the well-known Drude 
form: 



02£>(w) 



D 



TT 1 — ILOT 



(31) 



On the right-hand side of Eq. (|3lj) we introduced the 
(two-dimensional) Drude weight D = Y^u^u- ^ ne con ~ 
tributions of the partially filled subbands are 



D„ 



AttH 



2 Z v M v {e* v ) 
e tt- 



h' 2 

Here, e* is defined by the equation E v (e*) = and 
d 2 k 



(32) 



(2TTf 



\Ve k \'5(e-s k ) 



(33) 



is the weighted density of states4i The width 1/r of the 
Drude peak is obtained by the following weighted sum 



I = V— — 



(34) 



where t v is the relaxation time of the subband v. 

Experimentally, the Drude weight and the plasma fre- 
quency lj p d are usually determined by fitting the com- 
plex dielectric function to a Drude+Lorentz form^ In 
order to compare the calculated total Drude weight D = 
D v of the two-dimensional electron system with the 
plasma frequency uj p d measured in optical experiments 
on (MI) jv/(BI) a/ superlattices we assume independent 
quantum wells that yields the relation 



pD 



D N 
aN N + M' 



(35) 



(T2D = $(0) 



ATT 2 h 



S F L. 



(36) 



Sf = Y^u Sv(fy i s the total area of the Fermi surface and 
the mean free path is obtained as a weighted average, 



i= s^E^ (OK (0M0). 



(37) 



With the parameters estimated in the last section for 
the (LaTi03)jv/(SrTi03)M=io superlattices, we find an 
averaged mean-free path of L w 10 a which is about half 
of the value given in Ref. I20I . 

The electronic contribution to the thermal conductiv- 
ity follows from Eq. (|23|) . 



"2D 



(38) 



and is related to the electrical conductivity by the 
Wiedemann-Franz law typical for elastic scattering^ 



C. Thermopower or Seebeck coefficient 

Next we consider the thermopower Q, which, in an 
open circuit, is defined as the constant between electro- 
chemical field and temperature gradient E = QVT. For 
a metal at low temperatures, it follows from Eq. ([22l that 
Q is given by^ 



kir a , 



(39) 



E=0 



In the quantum-well system, we can interpret the ther- 
mopower as a weighted sum of the contributions of the 
partially filled quasiparticlc subbands, 



Q = 



J2u Qv&v 



a 



(40) 
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where cr„ and Q v are the conductivity and the ther- 
mopower of the subband v 1 respectively. Alternatively, 
one considers a weighted sum of the individual layers, 



Q = 



In the quasiparticle description, it is more natural to use 
the first interpretation and, in the following, we will focus 
on the multi-subband aspect. 

When a rigid-band picture is applicable, meaning that 
the dispersion E^ v is independent of the filling of the 
band, one can show that a multi-band system yields a 
small magnitude of Q. It is therefore less efficient com- 
pared to the single-band system with the best proper- 
ties because of possible cancellations in the enumerator 
of Eq. (j4"0)) due to electron-like (Q v < 0) and hole-like 
(<2„ > 0) contributions [see, e.g., Ref. However, this 
conclusion breaks down in strongly-correlated electron 
systems because the rigid-band picture is in general not 
applicable and therefore we believe that a multi-band sys- 
tem is not a priori less efficient. In fact, the effect of onsite 
correlations is strongest near half filling which leads to 
several overlapping subbands in superlattice structures. 



1. Thermopower: subband contributions 

The value of the thermopower for an individual sub- 
band is 



k B k B T 
e Z v 



r v {el) M v {et) 



z v 



(42) 



where r v {e) = t v (E v {e)) and the prime denotes the 
derivative with respect to e. There is an overall reduc- 
tion of the energy scale by Z v . As a result, the low- 
temperature slope of Q u is enhanced by a factor of 1 jZ v 
due to correlation effects J£ In the following we will dis- 
cuss the different contributions to Q v . 

The first term inside the square brackets of Eq. (|4"2"|) 
describes the influence of the scattering process. Usually, 
when considering elastic impurity scattering, this term is 
small and can be neglected. The second term describes 
the contribution from the uncorrclatcd band structure. 
It is sizable if the subband occupation n v is small, be- 
cause the particle-hole asymmetry is large near the sub- 
band edge. In fact, for the almost empty subbands with 
n u a 2 <C 1, we find 



Af v (et) ~ U + , 



1 



(43) 



On the other hand, in the presence of strong electron 
correlations, the third term Z' V IZ V is dominant for most 
subbands because the spatial non-uniformity of the sys- 
tem leads to a sizable logarithmic derivative of the quasi- 
particle weight at the Fermi energy 17 ' 47 

The dominance of the third term becomes evident in 
a study of the thermopower as onsite correlation effects 



(41) 1 
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FIG. 8: The thermopower Q„ and the relative conductivity 
ov/cr of the individual subbands as function of U r /U^ H for 
the N = 8 quantum well. The numbers on the right-hand 
side indicate the subband index in descending order in the 
subband occupation starting with v = 1 with the highest 
subband occupation. 



are varied. Q v is shown in the left panel of Fig. ([8]) as 
function of U r /Uf K for the N = 8 heterostructure at 
ksT = 0.01£. For the strongly renormalized subbands 
[v = 1 — 9 in Fig. [5]), we observe a clear qualitative 
difference between Hartree (U r < Uf R ) and Mott (U r > 
U^ R ) regimes. In particular, in the Mott regime, there is 
a noticeable spread of Q v mainly due to the third term 
in Eq. (|42"]) . On the other hand, the thermopower of the 
weakly renormalized subbands (y = 10 — 13 in Fig. [5]) 
is governed by Eq. (|4"3"|) and the change as function of 
U r /Uf K is due to the change of the subband occupation 
n v . Since t^/tv <ti 1, the qualitative behavior of the 
subband thermopower Q v is almost independent of the 
value of the impurity strength Vq. 



2. Dependence of the total thermopower on the impurity 
potential 

It is crucial for the total thermopower how the differ- 
ent contributions Q v are weighted. According to Eq. ([40]) 
the weighting factor is cr^/a, which depends on the Drude 
weight D v and on the relaxation time t„ . As an example, 
(7 j, I a is shown in the right panel of Fig. fusing the Born 
approximation of Eq. (|30j) corresponding to Vq — * 0. In 
this case, the most dominant contributions to Q in the 
Mott regime arise from the highest partially filled sub- 
bands (y = 10 — 13 in Fig. [5]). The reason is twofold: first, 
electron-like and hole-like contributions of the correlated 
subbands partially cancel, and second, the relative con- 
ductivity a v /<j of these bands is significantly reduced. 

However, for larger values of Vb, the contributions of 
the subbands subject to strong correlation effects become 
more and more important. Figure [5] shows the depen- 
dence of Q on the potential Vq for N = 5. For Vq — ► 
the weakly renormalized subbands are dominant, leading 
to a relatively large negative thermopower. For Vq — » oo 
the thermopower increases because the positive contribu- 
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FIG. 9: The thermopower Q of the N — 5 heterostructure as 
function of the impurity potential Vb, which governs the rela- 
tive contributions of the individual subbands. Also indicated 
is the value of Q when assuming a constant relaxation time. 



tions of some of the correlated bands gain more weight 
(see also Fig. [5]). Interestingly, for an even N, we observe 
that Q can take positive values at low temperatures be- 
cause there is one subband close to half filling with a 
dominant positive contribution. As a consequence, Q os- 
cillates as function of N (see Fig- EJ) - Moreover, the larger 
the N the larger the effect of cancellation becomes, which 
reduces the amplitude of the oscillation for N — > oo. 
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FIG. 10: The thermopower PF= aQ 2 as function of the in- 
verse width 1/N of the Mott-insulating material at fcsT = 
O.lt. The chosen parameters are as follows: the onsite inter- 
action is U r = 22t, the width of the band-insulator is M = 10, 
and the impurity scattering in the first Born approximation 
is parameterized by (Vb/t) 2 Ci mp =0.2. 



of Q 2 by increasing N (see Fig. [3]) and the volume term 
~ 1/(N + M), which reduces the total conductivity a. 



VII. CONCLUSIONS 



D. Thermoelectric application 

For thermoelectric applications, the dimensionless fig- 
ure of merit ZtT = aQ 2 T/{n e + n L ) should be as large 
as possible. This quantity provides a measure of the ef- 
ficiency of a material used for cooling or heating, respec- 
tively, as a thermoelectric power generator. It involves 
the static transport coefficients of the system, the elec- 
trical conductivity a, the thermoelectric power or See- 
beck coefficient Q and the total thermal conductivity 
k = K e + k l . A high value of a is necessary to minimize 
Joule heating, while a low value of K helps to maintain 
a large temperature gradient. In metal-oxide based ther- 
moelectric materials, the phonon thermal conductivity 
usually plays a dominant role and often n L 3> n e (see, 
e.g., Ref . [48l and l28h . Therefore, we will focus on a quan- 
tity involving only the electronic transport properties, 
namely the product PF= aQ 2 called the power factor. 
It is a common approach to optimize PF by varying ex- 
ternal parameters, and we will discuss the possibility to 
optimize PF by changing the superlattice parameters A^ 
and 

Let us, for example, fix M = 10. Since both Q and 
a depend also on N one can try to optimize the power 
factor by varying N. For the case of weak impurities this 
is shown in Fig. [10] for a temperature ksT = O.lt. We 
find that for this temperature the power factor takes the 
largest value at A" ~ 10. The maximum, as function of N, 
arises due to the competition between an enhancement 



In this work we have focused on the metallic properties 
observed in Mott-insulator/ band- insulator heterostruc- 
tures. We propose an analysis based on a quasiparti- 
cle picture that is expected to be valid in the Fermi 
liquid regime below the coherence temperature. A self- 
consistent renormalization of the quasiparticles was ob- 
tained by applying the slave-boson mean-field approx- 
imation and the quasiparticle transport was calculated 
in the relaxation-time approximation of the Boltzmann 
transport equation. 

The calculated optical conductivity is in good agree- 
ment with previously published spectroscopic measure- 
ments. Furthermore, the comparison with experiment al- 
lows to extract some of the key parameters of the model. 

We have shown that the correlation effects enhance the 
thermoelectric response as compared to a metal/ band- 
insulator interface. Furthermore, we have demonstrated 
that the design of the superlattice offers in principle ex- 
perimentally accessible parameters to optimize the power 
factor for a given temperature. However, the thermoelec- 
tric response, in particular the Secbeck coefficient, is very 
sensitive to the considered scattering mechanism. Even 
when focusing on s-wave impurity scattering at low tem- 
peratures, two completely different low-temperature be- 
havior of Q can be obtained in the Born and the unitary 
limit. Especially, this difference is manifest when consid- 
ering Q as function of the width N of the quantum well. 
On the other hand, measurements of the thermopower 
would actually be a good probe to test to which extent 
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the quasiparticlc description holds and, where applica- 
ble, to gain information on the delicate point of which 
subband contributions are most dominant for the quasi- 
particle transport. 

Whether the thermoelectric properties of a superlattice 
of the kind considered in this paper can overtop those of 
a bulk solid solution with the same composite remains 
unclear and further studies are necessary to clarify this 
point. In particular, we have neglected aspects concerned 
with the orbital and spin degrees of freedom, which may 
lead to ordered phases modifying the simple picture given 
here. Furthermore, our results are limited to the Fermi- 
liquid regime with coherent quasiparticlc transport. It 
is possible that the qualitative behavior at higher tem- 
peratures is quite different. For example, it is known 
that the thermopower in the single-band Hubbard model 
may even change sign with increasing temperature indi- 
cating the transition from electron-like coherent to hole- 
like incoherent transport. 22 ! 50 Similar effects may also oc- 
cur in the model discussed in this paper. Furthermore, 
for higher temperatures finite-lifetime aspects due to the 
electron-electron interactions and due to the coupling to 
the lattice degrees of freedom play a cumulative role. In 
this case a more sophisticated analysis is necessary. 
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APPENDIX A 

In order to clarify the approximations undertaken in 
the Boltzmann description for the dc transport coeffi- 
cients, we derive the transport distribution function &(E) 
starting from the linear-response Kubo formula^ and as- 
suming the validity of the quasiparticlc description out- 
lined in Sec. IIII Al This requires the evaluation of the 
current-current correlation function 



ir(ioj) = — 



T T j(r)-j(0) 



(Al) 



where Q = | a 2 is the sheet volume. A subsequent 
analytical continuation icu — > u> + i5 yields the sheet dc 
conductivity by taking the following limit 

"Im7r R (w)" 



(J2D = - lim 



(A2) 



The remaining transport coefficients may be obtained by 
applying the Jonson-Mahan theorem^ which holds for 
the model [Eq. |T|)] when replacing the long-range in- 
teraction part by a (self-consistently determined) one- 
particle potential^ 2 - The current operator for this model 



and for transport in the x,y direction is (h = 1) 

3 = - e X] V£kC ^ C ik<T' 



(A3) 



k is a two-dimensional crystal momentum and / the layer 
index. In the following we show that when neglecting ver- 
tex corrections in Eq. (|Alj) . one arrives at the transport 
distribution function [Eq. (|16[) ]. Thus, we approximate 



ir(iu}) 



e 



'I'll 



where the trace has to be performed over the different lay- 
ers and Q is the full Matsubara Green's function including 
all self-energy corrections. After letting iu> — > u> + i<5, one 
finds the form given in Eq. (|2~lj) . 



&2D 



dE 



1 1 < 



with the transport distribution function 



HE) 



ire 



£iv e 



k| 2 Tr 



i(k,£) : 



(A4) 



(A5) 



Here, we have introduced the spectral density matrix 
A(k,E) = — -Im G R (k, E) for which we assume (sec 
Sec. IIII Ap the following low-energy form 



A(k,E) 



^ (£-£k„) 2 +7L 



(A6) 



The quasiparticle life-time is given by the relation 71^ 
l/2T„(.Eki/), where 



7kjy = 



I 



(A7) 



In lowest order in temperature we can restrict the analy- 
sis to the Fermi surface. In this case 71^ is proportional 
to the impurity concentration Ci mp for dilute impurities. 



For non-degenerate subbands and in the limit ci, 







only terms, which are diagonal in the subband index, con- 
tribute to the trace in Eq. (|A5|) . After using the relation 
[Eq. p4p ] for the renormalization amplitude one ar- 
rives at Eq. (jT5J) obtained from the linearized Boltzmann 
equation (for to — 0). In actual calculations we estimate 
7k;/ by the approximations discussed in Sec. IIVBI 

In summary, we assume the validity of a Fermi-liquid 
description with local self-energy corrections, dominant 
s-wave scattering on dilute impurities, and that effects of 
weak localization can be neglected. Under these condi- 
tions we expect that, referring to the homogenous and 
isotropic system^ the transport distribution function 
used in the present work captures the main features in 
lowest order in temperature. 
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